JLAB-THY-03-199 
WM-03-110 

Higher moments of nucleon spin structure functions in heavy 
baryon chiral perturbation theory and in a resonance model 

Chung- Wen Kao 

Department of Physics and Astronomy, 
University of Manchester, Manchester, M13 9PL UK * 

^ . D. Drechsel, S. Kamalov 

O 



O 

(N 



Institut fur Kernphysik, Johannes Gutenberg- Universitdt, D-55099 Mainz, Germany 

M. Vanderhaeghen 



O 

Q 

qq ! Department of Physics, College of William and Mary, 



Williamsburg, VA 23187, USA and 

> 

CN . Jefferson Laboratory, Newport News, VA 23606, USA 



(Dated: 11th March 2009) 



cn ■ Abstract 
O. 

The third moment d<i of the twist-3 part of the nucleon spin structure function g2 is generalized 



to arbitrary momentum transfer Q and is evaluated in heavy baryon chiral perturbation theory 



(HBChPT) up to order 0(p 4 ) and in a unitary isobar model (MAID). We show how to link d<i as well 
as higher moments of the nucleon spin structure functions g\ and g2 to nucleon spin polarizabilities. 
We compare our results with the most recent experimental data, and find a good description of 
these available data within the unitary isobar model. We proceed to extract the twist-4 matrix 
element /2 which appears in the 1 /Q 2 suppressed term in the twist expansion of the spin structure 
function g\ for proton and neutron. 
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I. INTRODUCTION 



The forward scattering of spacelike virtual photons (with virtuality Q 2 ) on the nucleon, 
allows one to study sum rules which relate nucleon structure quantities to inclusive elec- 
troproduction cross sections. At large Q 2 , it yields the sum rules studied in deep- inelastic 
scattering (DIS) experiments (see Ref. |]J for a recent review). The study of such sum rules 
as function of Q 2 from the real photon point to large Q 2 , opens up the perspective to in- 
terpolate between the non-perturbative and perturbative regimes of QCD, as one goes from 
low to high Q 2 . 

The moment d 2 of the nucleon spin structure functions can be measured by scattering 
longitudinally polarized electron beams off nucleon targets with transverse and longitudinal 
polarizations. Being a higher moment in the Bjorken variable x, d 2 contains appreciable 
contributions from the resonance region. It is therefore the aim of this letter to study the 
threshold and resonance contributions to d 2 within the frameworks of heavy baryon chiral 
perturbation theory (HBChPT) and of a unitary isobar model (MAID). 

Since a measurement of d 2 requires also transverse polarization, experimental information 
on this observable has become available only recently at SLAC P, [J and JLab [4] . Further 
experiments are underway or proposed at JLab [Mez03]. In particular, the SLAC experi- 
ments |3| yielded the values : d\ = 0.0032±0.0017 and d% = 0.0079±0.0048 at Q 2 = 5 GeV 2 . 

At large momentum transfers, a non-zero value for d 2 directly measures a twist-3 quark- 
gluon matrix element in the nucleon, which has been estimated within several model cal- 
culations as well as within lattice QCD. In this paper, we show that when generalizing the 
definition of d 2 to arbitrary momentum transfers, its physical interpretation can be expressed 
in terms of generalized (i.e. Q 2 dependent) spin polarizabilities of the nucleon. These gen- 
eralized polarizabilities appear in the analysis of the forward virtual Compton scattering 
amplitude (for a review, see Ref. M). Using the calculation of these spin polarizabilities 
within HBChPT j^J (see also Refs. £j, 0]), and within a unitary isobar model we eval- 
uate d 2 . In an analogous way, the third moment of the spin structure function g\ can be 
expressed through these spin polarizabilities. These generalized definitions for d 2 and the 
third moment of g\ constitute useful observables to interpolate between a hadronic descrip- 
tion at low Q 2 and a partonic description, based on the Operator Product Expansion (OPE) 
at large Q 2 . 
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Besides the twist-3 matrix element d 2) the l/Q 2 suppressed term in the twist expansion 
of the first moment of g\ also contains the twist-4 matrix element f 2 . It has been suggested 
(see e.g. Ref. Q]) that the matrix elements d 2 and f 2 enter in the response of the color 
electric and magnetic fields to the polarization of the nucleon in its rest frame, which can 
be expressed in terms of gluon field polarizabilities. Therefore an extraction of f 2 , besides 
d 2 , can yield us interesting new nucleon structure information. Such an phenomenological 
extraction was peformed for the first time in Ref. |h| based on the SLAC data of Ref. Q]. 
With the advent of the recent JLab data for the first moment of the proton and neutron spin 
structure function gi in the intermediate Q 2 range, we may now perform a phenomenological 
analysis based on all available data to extract the twist-4 matrix element f 2 . 

The outline of the paper is as follows. In Section |H1 we review the twist expansion of 
the moments of the nucleon spin structure functions g\ and g 2 and show how the definitions 
for the third moments of the nucleon structure functions g% and g 2 can be generalized at 
arbitrary Q 2 in terms of nucleon spin polarizabilities. 

In Section we use the calculation of these spin polarizabilities within HBChPT to 
evaluate d 2 and the third moments of g\ and g 2 . 

In Section IIV1 we show how these observables can be expressed in multipoles which arise 
naturally when performing an evaluation within an isobar model. 

In Section we show our results for d 2 and the third moments of g\ and g 2 within both 
ChPT and a unitary isobar model and compare with recent experimental data for d 2 . We 
also compare our results for the first moment of g\ with the available data for both proton 
and neutron and use the twist-expansion for the first moment of g\ to extract the magnitude 
of the twist-4 matrix element f 2 . 

Finally, we give our conclusions in Section IVT1 



II. TWIST EXPANSION AND MOMENTS OF NUCLEON SPIN STRUCTURE 
FUNCTIONS 

At sufficiently large momentum transfer Q 2 » Aq CD , one can perform a twist expansion 
for the moments of the nucleon structure functions. We will consider in this work the 
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moments of the nucleon spin structure functions g\ and g 2 , which are defined as : 

rf } (Q 2 ) = f 1 dxx n ~ 1 g 1 (x, Q 2 ) , n=l,3,5,... (1) 

Jo 

ifV) = f dxx n - 1 g 2 (x, Q 2 ), n=l,3,5,... (2) 

JO 



In particular, for the first moment of g\ such a twist expansion can be written as |I2|,|14| : 

r?(Q 2 ) = rfUW 2 ) + ^ (a 2 (g 2 ) + 4d 2 (Q 2 ) + 4/ 2 (g 2 )) + oQ) . (3) 

The remaining Q 2 dependence in the coefficients Ti\ w _ 2 , a 2} d 2 and f 2 of the twist expan- 
sion is logarithmic. In particular, the three-loop result for the leading (i.e. twist-2) term 
^x\ w - 2 (Q 2 ) is given, for 3 quark flavors, as 

f[\w-2{Q 2 ) 
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(4) 



where in the term proportional to g& the sign + (-) corresponds with proton (neutron) 
respectively, g& is known from neutron beta decay, and a§ is extracted from hvperon weak 
decay data assuming SU (3) symmetry. The updated values for g^ [3] and as 3| are given 
by : 



g A = 1.267 ±0.003, 
a 8 = 0.585 ±0.023. 



(5) 
(6) 



Furthermore, in Eq. (J3J), Oo° = clo{Q 2 = oo) is the flavor singlet axial charge, which is fixed 
by measuring the first moment 1^ at a sufficiently large scale. 

The 1/Q 2 suppressed terms in Eq. (JHJ) are of three different natures. The term propor- 
tional to a 2 {Q 2 ) arises due to target mass corrections and is given by the twist-2 part of the 
third moment of gi : 



a 2 (Q 2 ) = 2Tf\ w _ 2 {Q 2 ). 



(7) 



The terms proportional to d 2 and f 2 in Eq. Q correspond with dynamical higher-twist 
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corrections. The function d 2 corresponds with the matrix element of a twist-3 quark gluon 
operator, and can be expressed in terms of the twist-3 part of the spin structure function g 2 



as 



d 2 (Q 2 ) = 3 £ dxx 2 [g 2 (x, Q 2 ) - gf w {x, Q' 



(8) 



n 



where g^ w is the twist- 2 (Wandzura- Wilczek) part of g 2 [20] . Using the Wandzura-Wilczek 
relation, one can express Eq. (jHJ) as : 



d 2 (Q 2 ) = J q dxx 2 (3 g 2 (x, Q 2 ) + 2 gi (x, Q' 



(9) 



Several model estimates as well as lattice QCD calculations 



22j have been performed for the 



twist-3 matrix element d 2 . In particular, an estimate in the instanton vacuum approach 



2i|, 



where d 2 is parametrically suppressed due to the diluteness of the instanton medium, pre 



dieted d 2 to be of order 10~ 3 . Also a revised lattice calculation [22j supports this small 
value. Such a small value for d 2 , of order 10~ 3 , is in agreement with the recent experimental 
results from SLAC 0|. 

The function f 2 in Eq. corresponds with the matrix element of a twist-4 quark gluon 
operator. Unlike a 2 and d 2 , the matrix element f 2 cannot be expressed directly in terms of 
moments of g\ and g 2 as in Eqs. (JJJ) and Q. It can however be extracted phenomenologically 
from the twist expansion of Eq. Q as : 



4M| 
9Q 2 



/ 2 (Q 2 



(i) 



Ml 
9Q 2 



a 2 (Q 2 ) + Ad 2 (Q 2 



(10) 



by using the experimental information on the full Q 2 dependence of and provided one 
knows the twist-2 and twist-3 contributions r[ tw ^. 2 i a 2 and d 2 . In this way, f 2 has been 
estimated in Refs. for Q 2 > 0.5 GeV 2 where the twist expansion was assumed to 

hold. Several model calculations have also been given for f 2) _e.g. in the framework of the 



21] and using QCD sum rules j^J, to which wo refer further 



instanton vacuum approach 
on. 

Although the Operator Product Expansion (OPE) of Eq. (jHJ) is only defined for Q 2 >> 
Aq CD , Eq. (jHJ) can be used to define d 2 outside this range by keeping the full Q 2 dependence 
in the structure functions which appear in the integrand. Analogously, we can study the 
Q 2 dependence of the higher moments of g\, such as r x \Q 2 ), outside the range of validity 
of the OPE. We next show that at low Q 2 , the physical meaning of d 2 and can be 



expressed in terms of nucleoli polarizabilities which will be calculated in this work within 
Chiral Perturbation Theory and estimated within a phenomenological resonance approach. 
At large Q 2 , d 2 and tend into the matrix elements appearing in the twist expansion of 
Eq. © and display only a logarithmic Q 2 dependence. Defined in this way, d 2 and 1^ are 
useful observables to interpolate between a hadronic description at low Q 2 , involving the 
polarizabilities of the system, and a partonic description based on the OPE at large Q 2 . 

Starting with d 2 , one can split the integral of Eq. © into an elastic contribution at x — 1 
and an inelastic contribution. The elastic contribution is given by : 

While this contribution vanishes like Q~ 8 for Q 2 — > oo, it increases with smaller values of 
Q 2 and approaches df (0) = Ge(0) ■ Gm(0). The inelastic contribution to d 2 corresponds 
with the integral over the excitation spectrum : 



d inel 



(Q 2 ) = J dxx 2 [3g 2 {x, Q 2 ) + 2 9l (x, Q 2 )) . (12) 



where Xq = Q 2 / (2M^m n + m 2 + Q 2 ) is the threshold for pion production. 

To evaluate Eq. (fT2*|) . we can equivalently express the third moment of the twist-3 part 
of g 2 in terms of the spin-dependent doubly virtual Compton scattering amplitude in the 
forward direction (VVCS). Following the notations of Refs. (| Q], we use the VVCS am- 
plitudes gTri^iQ 2 ) and gLT(v,Q 2 ), where v is the lab energy and Q 2 the virtuality of the 
virtual photon, and T (L) denotes the transverse (longitudinal) virtual photon polarization. 

The imaginary parts of gxr and g^T are related to the virtual photon absorption cross 
sections ott and olt-, multiplied by a photon flux factor K (with dimension of energy) [50]. 
These partial cross sections are related to the nucleon structure functions g± and g 2 as : 

K-a TT = ^^(^(x,Q 2 )- 7 2 ^(^g 2 )), (13) 

K-a LT = ^^l( gi (x,Q 2 )+g 2 (x,Q 2 )y (14) 

with 7 = Q/v and x = Q 2 /(2M N v). 

For the non-pole (i.e. inelastic) contributions to gxr and giT, one can perform a low 
energy expansion (LEX) as follows 



Re g TT{ y : Q2) _ Re gfffa Q2 ) = /A( g 2) v + 7o(g2) u 3 + 0(l/5) j (15) 



Reg LT (v,Q 2 )-Regi¥(v,Q 2 ) = (-^\qi 3 (Q")+Q5 lt (Q 2 )^ + 0(^).(1Q) 



For the 0(y) term in Eq. (|15j) . one obtains a generalization of the GDH sum rule as : 

4 7r 2 a em Juq v v 



2 M% r*o 



Q2 jf° <fc |^ Q 2 ) - M x * ^ (x? Q2) I ) (17) 

and recovers the GDH sum rule at Q 2 = 0, as Ia{0) = ~ K /v/4, with the nucleon anoma- 
lous magnetic moment (n p = 1.79, n n = —1.91). 

Furthermore we introduce the integral I\, which is related to the inelastic part of the first 



hq 2 ) = ^ r dx 9l (x,Q 2 ), (is) 



moment of g\ as 



Q 2 Jo 

and which also approaches the GDH value at Q 2 = 0. 

The higher order terms in Eqs. (fT3j) and (fTBj) can be expressed in terms of nucleon spin 
polarizabilities, see Ref. p. In particular, the 0(u 2 ) term in g^ is given by : 

= aem ^ 6 6 M|r £° d» 2 {^i (x, g 2 ) + s 2 (x, g 2 )} . (19) 

Combining Eqs. (|17I18I19)) . we find that the inelastic contribution to the third moment 
g?2 of Eq. (|T2J) can be expressed as : 

™ - m {-ah? (w 2 > - ^ 2 >) + i fe « 2 >} ■ < 20 > 

We can perform a similar analysis for the moments 1^ and defined through Eqs. (JIJ 
and (J2J). The elastic contribution to 1^ is given by : 



While this contribution vanishes like Q~ 8 for Q 2 —>■ oo, it approaches 1^ 6 '(0) = 1/2 ■ 
Ge(0) ■ Gm(0) at the real photon point. The inelastic contribution to corresponds with 
the integral over the excitation spectrum : 

V (n)inel {Q2) = p dx x n-l Q 2 } (22) 



The third moment Y^ mel can be expressed in terms of the quantities introduced in 
Eqs. (HHHm as : 



if *-(«») - ^{m^('^-'^) + dr MQ2) > ■ (23) 



For T 2 , the elastic contribution is given by : 

rt )el (Q 2 ) = -\ A J\ ^ (Gm(Q 2 ) - Ge{Q 2 )) Gm{Q 2 )- (24) 



2 AM 2 N + Q 2 

This contribution vanishes at the real photon point. The inelastic contribution to 
corresponds with the integral over the excitation spectrum : 



-p(n) inel ( f)2 
1 2 



(Q 2 ) = / dxx n - 1 g 2 (x, Q 2 ). (25) 
Jo 



The third moment Y 2 3 ^ mel can be expressed in terms of the quantities introduced in 
Eqs. (IT7lrH|) as : 

r (3W (g2) = ^ (g2) _ Ji(Q2) j (26) 



Note that is linearly dependent on d 2 and as it can be expressed as : 

rf } (Q 2 ) = l(d 2 (Q 2 )-2Tf\Q 2 )) . (27) 

Although we restrict our investigation in this work up to the third moment of g\ and g 2 , 
one can in principle extend the analysis to the higher moments and r 2 , with n = 5,7.... 
Their expressions involve hig her spin polarizabilities, as have e.g. been introduced at the 
real photon point in Ref. |25j |. 



III. CALCULATION OF THE MOMENTS d 2 , if \ AND r| 3) IN HBChPT 



Eqs. (J! 



show that d l 2 nel {Q 2 ), mel (Q 2 ) and T ( 2 >meL (Q 2 ) can be expressed in terms 



-,(3) inel . 



of Ia, h and 5lt appearing in a low energy expansion of the forward VVCS amplitudes . 
These expressions have been calculated at low Q 2 in HBChPT, which allow us to construct 



the HBChPT predictions for d™ 1 , rS 3)me ' and T ( 2 )inel 
HBChPT have been obtained up to 0(p 4 ) in Refs. 

1 
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The expressions for I a and l\ in 



h(Q 2 ) 



16 



[(k s + K V T 3 ) 



+ 



9A M N m n 

(4ttF^) 2 ' M N ' 32 



{(_10 - 12k,) + (-2 - 12«.)t 3 

1 uu 

+ [(20 + 24ac) + (4 + 24k s )t 3 ] ■ - tan _1 [-] 

7/1 

+ [(3 + Qk v ) + (3 + 10k s )t 3 ] • lutan-^-]} , 



Ia(Q 2 ) = —[(«. + k.ts) 3 ] 
lb 



+ 



9a^n o /v^+4 



(4tt^)2 



sinh- 1 !-] - ll 



+ [(20 + Ak v ) + (4 - 12k s )t 3 ] ■ - tan" 1 [ 



w 



+ [(3 + 3k v ) + (3 - /? s )r 3 ] • wt&n- l [-)} 



(28) 



(29) 



with iy = y S". These expressions are given in terms of the renormalized isoscalar (k s ) and 
isovector (n v ) anomalous magnetic moments, whose physical values are given by k s = —0.12 
and k v = 3.70. Furthermore, throughout this paper we use the values qa = 1.267, F n = 
0.0924 GeV, and m v = 0.14 GeV. While both I x of Eq. and I A of Eq. flU approach 
the GDH value of —k%/ at Q 2 = 0, their slopes at this point differe substantially : //(O) = 
(6.31 + 0.66 r 3 ) / GeV 2 and 1^(0) = - (12.43 + 2.09 r 3 ) / GeV 2 . 

The longitudinal-transverse generalized forward spin polarizability Slt of Eq. (fTKJ) has 
also been calculated in HBChPT Q (see Ref. Q for the corresponding calculation within 
the framework of relativistic baryon ChPT). For the 0(p 3 ) term, the result is : 



LT 



(Q 2 



2 

®-em9A 



mi 



. sinh 1 [— 1 

3wVw 2 + 4 L 2 J 



(30) 



and the 0(p 4 ) correction is given by : 



'LT 



(Q 2 



2 

QemQ A 



N 



7T 

64 



{(-16 + 8n,) + (-8 + 16k s )t 3 
+[(-M + 8k v ) + (-6 + 8k s )t 3 ] 
+ [(-9 - 12/c) + (-9 - 4« a )r 3 ] 



1 

w 2 

1 , 
— tan I — 

w L 2 



w , 



w , 



+ [(-84 - 16k,) + (12 - 16«.)r 3 ] ■ tan"^- 



w 3 1 2 



■[4-(12 + 16« a )r 3 ] 



4 + w 2 



+ 128 



3 + w 2 
Aw 2 + w A 



X31) 



9 



Using Eqs. (|2*8*l - IBT ]) . we can now construct the result of d™ el from Eq. (j^UJ). The result 
at 0{p 3 ) for df el is 



(Q 2 



9l 



8M 2 N (AnF^m 2 V3w 3 V^T4 
and the 0(p 4 ) correction to d™ el is given by : 



( -24 -2w 2 2 
I — . , . smh x [-] + — 



(32) 



9l 



" 8M 2 N (AnF^ml M N 32 
x {[(-8 + 4^) + (-4 + 8« s )t 3 ] 



+ 



;-57 - 12k v + 16w 2 ) + (-9 + 4k s )t 3 



w- 



+ - [(9 + 12^) + (9 + 12rc s )r 3 ] • - tan" 1 & 

— — 



+ [(18 + 24/^) + (18 - 8k s )t 3 ] • -3 

3 + w 2 , 
2 + 64 (6 + 8k s )t 3 



tan 



w- 



A + w 2 



(33) 



Analogously, we can also construct the results for Y^ mel and T 



t(3) inel 



of Eqs. 



in 



HBChPT by using the expressions of Eqs. (1281 -13 1 )1 . The result at 0(p 3 ) are given by : 



(S)O(p 3 ) 



,(3)0(p 3 ) 



(Q 5 



(Q 2 ) 



gl ( 12 + 4 " 2 shlll 

g 6 <? 2 / 4 + w 2 



-i 



AM 2 N {47iF n ) 2 m 2 n Vw 3 V^T4 
and the 0(p 4 ) corrections are given by : 



2 



w- 



sinh *[-] - 



(34) 
(35) 



(S)O(P 4 ) 



(Q 2 ) 



si 



7T 



8M^ (47rF^) 2 m 2 M N 64 
x {[(-40 + 4/<) + (-4 + 8<)t 3 ] 

+ [(33 + 36k,) + (9 + Ak s )t 3 ] ■ — 



+ - [(-45 - 60^ 



1 



UK 



(-45 - 36k s )t 3 ] ■ -tan- 1 [ 2 

I TV 

+ [(-162 - 72k v ) + (-18 - 8k s )t 3 ] ■ — tan~ x [ 



+ 



3 + w 2 . 
2 + 64 (6 + 8k s )t 3 



w- 



1 



A + w 2 



(36) 



F (3)0(P 4 ) 
1 2 



(Q 2 ) 



<7i 



7T 



8M 2 , (47rF T ) 2 m 2 Mjy 32 
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X 



-30 - 16k,, + 16w* 



6r 3 



w- 



1 i w 

+ [(9 + 12k„ + (9 + 8k s r 3 • - tan" 1 - 

w 2 

+ [(60 + 32^) + 12r 3 ]-^tan- 1 [^]) . 



(37) 



From Eqs. (|3*2l - IHTj) . one can extract that the following prediction of HBChPT in the limit 
Q 2 ->0: 

Q 6 1 



df el (Q 2 



%M 2 N 2a em 



;-7o(0) + 35 LT (0)) 



Q 6 



F? )inel (Q 2 



4&M 2 N {±>itF„) 2 ml 
Q 6 1 



8M 2 N 2a em 



7o(0) 



-,(3) inel 



(Q 2 



ASM 2 N (AnF n ) 2 ml 

Q 6 l 



Tfl 7T 

1 + TT ■ « t( 21 + 9k -) + (~ 6 + 14 ^) r 3] ^ • ( 38 ) 



{2 + ^ • I [(-15 - 3k„) + (-6 - k s )t 3 ]} , (39) 



SM 2 N 2 a em 

Q 6 



M N 

7o(0) + 5 Lr (0)) 



A8M 2 (47rF n ) 2 m 2 



Tfl 7T 1 

-1 + j^' -[(27 + 7k u ) + (6 + 6k s )t 3 ]| , (40) 



where 7o(0) is the forward spin polarizability at the real photon point, which has also been 
calculated in HBChPT to 0(p 4 ), see Ref. 0. 

We have also studied the contribution of the A(1232) resonance in the "small scale 
expansion" 0(e 3 ), which uses the quantity A = M^ — M^ as a further expansion parameter. 
Our results for the quantities Ji, I a and 5lt are : 
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siAQ 2 
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18 v A 
8Q 2 M 2 g 



ttAN 
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(4tt^)2 . 


-32a e „ 
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1 j 7v AN 
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(4ttF,)2 


16a em 


Aan 1 



9 (4ttF w ) 



1 , x 2 (l-2x) / 2 
r/.r — i — - 



A*o- 



ln\ji + Jl4 ~ 1] 



/^o - 1 



},(41) 



.2; 



cfe ^ '" 2 



K - 1: 







1 , x 2 (l -2x), o 

GSX n 



3/i, 



ln[/i + J/4- 1] 



^-1 



/o 







ln[// + a/Vo - 1] 



//§ - 1 



(42) 



with m = \jrn1 + x(l — x)Q 2 and /io = Furthermore, Ci and g nl \N are the leading 



order jAN and tcAN coupling constants respectively. In the large Nq limit of QCD, they 
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are related with k v and as : 

Gi 



2V2 



2V2 



9A- 



(43) 



Furthermore, the 0(e 3 ) A contribution to g?2 is 

Q 6 



N 



r 1 Gi 2 
\1S K M N ) 
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A 2 


32 9ttan 


/' 


27 (47TF,) 2 


Vo 


, 24 9wAN 


/' 


9 (47TF.) 2 


Jo 



X 



dx — ^(/Ug — 1) 



-2 



rn, 



ln[/i + J/4 ~ 1] 



0o - 1 



r/x X2(1 . 2X) (^ - I)" 1 



1 - 0o 



nl - 1 



(44) 



and the 0(e 3 ) A contribution to Y^ mel is 



rf )A (Q 2 ) 



Q ( 



-1, Gi 



8Af£ I 18 'M^v 



A 2 



16 <? 



27 (4^) 



* /■ 1 

tFJ 2 Jo mn 



-2 



o ln [0O + V0O-l] 

A*o + 2 — 3 0o- 

>o - 1 



(45) 



In the limit Q 2 — >■ 0, these quantities approache the values 



rf )A (Q 2 ) 



1 / K v \2 ^ 

8M^ ll6 l M7/ ' A 2 
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3 (47rF^) 2 m 2 
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2 (47rF^) 2 m 2 



8M^ I 16 V M 
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6 (47rF 7r ) 2 m 2 
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(0 2 - 1) 
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1-0 
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(0 2 " 1)- 



„ 2 + 2 - 3/, ln[/i + 



(46) 



(47) 



with n = Analogous formulae can be derived for I2 A (Q 2 ) by taking the combination 
of Eq. U23). 
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IV. MULTIPOLE CONTENT OF THE MOMENTS 



By use of Eqs. (I13I14|) . we can express the moment d<i by the partial cross sections ott 
and Of/r : 

g6 ™ irf,, r>2\ / q,,2 , or>2 



d™ el (Q 2 ) 



32ix 2 a em M 2 



em lvl N Ju 



v 2 (v 2 + Q 2 ) 



2 3j^+2Q^ 2 

-o-ttw Q ) + — ^ — vuriy, Q ) 



(48) 



and the multipole content follows by inserting the expansion of these cross sections in terms 
of the electric (E), magnetic (M), and Coulomb or scalar (S) multipoles 24 1 



Ka TT (lir) = x TT ^-(/ + l)[(/ + 2)(|^ + | 2 + |M i+1 ,_| 2 ) 
i 1 

-/(|£ m ,_| 2 + \M l+ \ 2 ) + 21(1 + 2)(E* l+ M l+ - £f +lj _M m ,_; 
= x TT (\E 0+ \ 2 + \M 1 .\ 2 -\M 1+ \ 2 + QE* 1+ M 1+ + 3\E 1+ \ 2 ± ...) , 



(49) 



Ka LT (lTi) 



I z 

• [s; + ((l + 2)E l+ + IM 1+ ) + S? +h _(lE l+1 ,„ - (I + 2)M m , 
xlt (So + E 0+ - S' 1 *_M 1 _ + 2S* + (M 1+ + 3E 1+ ) ± . 



where 



x TT = x TT (v, Q 2 ) = Air^J (u - u ) (u-u + 2m v 
xlt = %lt(v, Q 2 



\ 



1 + 2u/M N - Q 2 /M 



N 



(50) 



(51) 
(52) 



1 + v 2 /Q 2 

The multipoles are generally functions of the cm energy of the photon, Aii± (u>, Q 2 ), with 

M N u - Q 2 



2M N u ■ M 2 Q 2 



(53) 



In comparing with Ref. 



24j it should be noted that we have changed the sign of the 



partial cross sections in order to stay in accordance with the usual notation of deep inelastic 



scattering, i.e., (t^t — — dj'j' and &LT = —^'lti where a' TT and a' LT were used in Ref. j24j 
We also note that expressions like Sq + E q+ should be read as Re (Sq + E 0+ ). 

The expression for the moments , and /2 can be derived similarly as Eq. ()48j) . 

Q 6 , K(v ^ 2 

+ Q 2 ) 



Tf )me \Q 2 



327r 2 a Rm M 2 



em lvl N Jv 



/>ra(*MVi,M ! )) , (54) 
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^ {Q) = 32^M L ^ v^ + QZ) [~ aTT ^ Q ] + Q° LT ^ Q ] ) ' (55) 
finely 9Q 2 (1) 9Q 4 roo K(u,Q 2 ) f 2 Q 2 , 



4M 2 ' 32n 2 a em M 2 N J„ v 2 + Q 

Q 6 r, K(u,Q 2 ) ( 3(^ 2 + g 2 ) 2 \ 

647r 2 a em MXr ^0 z/ 2 (z/ 2 + Q 2 ) \ Qv 

V. RESULTS AND DISCUSSION 

In Fig. H we show the results for the Q 2 dependence of the moment d 2 . At large Q 2 , the 
values of c? 2 have been obtained for both proton and neutron by DIS experiments at SLAC 
using a transversely polarized target jjj. From these experiments, the following values have 
been extracted at an average value of Q 2 = 5 GeV 2 |j| : 

c^(Q 2 ^5GeV 2 ) = 0.0032 ±0.0017, 

d n 2 {Q 2 « 5 GeV 2 ) = 0.0079 ± 0.0048 . (57) 

In the low and intermediate Q 2 region, we show the results for c? 2 according to the generalized 
definition of Eq. (JHJ). It rises strongly at the lower Q 2 , displaying a Q 6 dependence at low 
Q 2 , and tending to a small constant value asymptotically, corresponding with the twist-3 
matrix element entering in the OPE of Eq. (J3J). Due to this structure, d 2 is an interesting 
observable to interpolate between a hadronic description at low Q 2 and a partonic description 
at large Q 2 values. One sees from Fig. ^ that the ChPT results rise strongly with Q 2 , 
and its Q 6 dependence at the low Q 2 values is given by Eq. (j38|) . The large difference 
between the 0(p 3 ) and 0{p i ) results originates from the known large difference between the 
0(p 3 ) and 0(j» 4 ) HBChPT results for the forward spin polarizability 70 [3, l^l l^of - It is 
interesting to note however that the 0(p 3 ) HBChPT result is in good agreement with the 
phenomenological MAID estimate up to Q 2 ~ 0.25 GeV 2 . The convergence of the heavy 
baryon chiral expansion for the forward spin polarizability has been recently investigated in 
Ref. P| using the Lorentz invariant formulation of baryon ChPT of Ref. j^|. It was found [l| 
that the main reason for the slow convergence of 70 in the 1/Mjv expansion is due to the 
slow convergence of the Born graphs. The corresponding one-loop relativistic calculation 
cures this deficiency. However the relativistic result to fourth order in the chiral expansion, 
even when supplemented with A and vector meson contributions, still does not agree with 
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the data, suggesting to envisage a fifth-order calculation for 70 in future work. 

Very recently, the inelastic contribution to d 2 has been measured for the neutron at 
i—a teC? ; valuesatJLaVHaHAQ.lt is seen from Fig. Q that the phenomenological 
resonance prediction using the MAID model for the neutron shows an excellent agreement 
with these data. We will discuss the MAID results for d% further on, to get some better 
understanding of the origin of this good description. 

In Figs. 121 and 01 we show the results for the Q 2 dependence of the third moments of the 
nucleon spin structure functions g\ and g 2 respectively. From Eqs. ()39I40|) . we see that the 
inelastic part to the moment Yx is proportional to Q 6 ■ 70 at small Q 2 , and that the inelastic 
part to the moment also contains 70 at small Q 2 . Therefore the HBChPT results for 
those moments directly reflect the poor convergence of the chiral expansion for 70 . The 
MAID model predicts a value for 70 at the real photon point as : 7o = —0.707 • 10 _4 fm 4 , 
which is in relative good agreement with the experimental value (for its determination, see 
Ref. Q) : 

7^ = [-1.01 ± 0.08 (stat) ± 0.10 (syst)] • 10~ 4 fm 4 . (58) 

Going to larger values of Q 2 , the MAID results for the proton change sign. At large Q 2 , the 
value of is known from DIS experiments as shown on Fig. |21 We also compare in Fig. |21 
the MAID result with the resonance contribution to 1^ (corresponding with W < 2 GeV) 
as extracted from the fit of g\ data. From this comparison, we see that the 7r-channel alone 
underestimates the total resonance contribution at larger Q 2 . For the proton channel, we are 
also able to provide an estimate for the rjN and ttttN channels which provide the dominant 
virtual photon absorption cross sections beyond nN. It is seen that the sum of the n + Tnr + r] 
channels nicely joins the W < 2 GeV part of the DIS parametrization (lower shaded band 
in Fig. |21) for Q 2 > 3 GeV 2 . To gain an understanding of the gradual transition in l j 3 \ 
from the resonance dominated to the partonic regime, we generalize a suggestion of Ref. 33j 
for the first moment I\ of g\ to parametrize its Q 2 dependence through a vector meson 
dominance type model. This model for Ii was refined in Refs. by adding explicit 

resonance contributions which are important at low Q 2 as seen above. We can generalize 
this phenomenological parametrization for the inelastic part to an arbitrary moment of gi 
through the interpolating formula : 

r (n)inel {Q2) = T (n)res {Q2) + ^ (Q*lf >" + f «.(»)) j (59) 
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where res (Q 2 ) is the resonance contribution evaluated using the MAID model, and /x is 
a vector meson mass scale (// ~ 0.77 GeV), which governs the transition to the asymptotic 
value as . The coefficient 

c (n) jg 

chosen in such a way that reaches its known value at 
the real photon point. For this is fixed by the non-resonant contribution to the GDH 
sum rule as : 

= ^ (-K 2 /4-/r(0)), (60) 

where 7[ es (0) is the resonance contribution (W < 2 GeV) at the real photon point to 1%. For 
the third moment, the coefficient is fixed by the very small non- resonant contribution 
to the forward spin polarizability : 

c (3) = l^S-^-(7b(0) -tHO)), (61) 
16 a em 

where 7o es (0) is the resonance contribution (corresponding with W < 2 GeV). Analogously, 
one can express the coefficients c- n > for higher moments (n = 5, 7, ...) through the very small 
non-resonant contributions to the higher spin polarizabilities of the nucleon. As these higher 
polarizabilities are practically completely dominated by the low-energy excitation region (i.e. 
resonance region), it is a very good approximation to take ~ for n > 5. 

In Fig. El we also show the result of the interpolating model of Eq. for T\ of the 
proton, by using the MAID model for the resonance contribution rj eJ and using as asymp- 
totic contribution for the proton : ~ 0.012. It will be interesting to compare these 
predictions with experimental results which can be extracted from corresponding analyses 
for the lowest moments of the proton [36] and neutron 3tJ structure function g\. 

Besides the inelastic contributions, we also show in Figs. 11121 for comparison the elastic 
contributions to d,2 and according to Eqs. I)ll|21j) . The elastic contribution is calculated 
by using for the proton the form 



data for the ratio of G v E jG 



actor parametrization of Ref. |38| including the new JLab 
For the neutron, the elastic contribution is calculated 



using for G \ the recent fit of Ref. (41] (following the Galster form), and for G\ the recent 
fit of Ref. |42j. We see from Fig. [I] that for the proton the elastic contribution largely domi- 
nates the inelastic one for both c? 2 and when Q 2 < 1 GeV 2 . For the neutron, the elastic 

(3) 

contributions to d,2 and 1^ are of comparable size to the inelastic ones and vanish at the 
real photon point. 

In Fig. HI we show the separate contributions of ott and <jlt to and r^ n to gain 
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some insight in the nature of the photon absorption mechanism involved. In the case of 
the twist-3 moment d 2 of Eq. (|48jh <jlt is enhanced by the kinematical term in front, and 
therefore becomes dominant at large values of v and Q 2 . The situation is quite different for 
the third moment n of Eq. (J54)) . in which case the LT term vanishes at Q 2 = and is 
suppressed for large v. 

The multipole content of d 2 is displayed in Fig. |3] for both proton and neutron. This 
moment is dominated by an interplay of s- and p-wave multipoles up to Q 2 ~ 4 GeV 2 . The 
striking difference between the proton and neutron is essentially due to the excitation of the 
Sii(1535) resonance. According to the Particle Data Group [Tj|, the helicity amplitudes for 
the 5n(1535) are A p 1/2 = (0.090 ± 0.030) GeV" 1 / 2 and A% /2 = (-0.046 ± 0.027) GeV" 1 / 2 , 
with the result that the cross sections for Su excitation on the proton are about a factor 4 
larger as c omp ared to the neutron. It is also known that this resonance has a very hard form 
factor 43j,|44|, and therefore this effect persists up to very high momentum transfers. This 
leads to the difference that d 2 changes sign around Q 2 ~ 1.25 GeV 2 , whereas d 2 remains at 
the same sign as confirmed by the data. 

Fig. El shows the multipole analysis for the moment T^ n . The conclusions are similar 
as in the previous case. The strong coupling of the Su to the proton leads to the distinct 
feature that changes sign at Q 2 ~ 1.25 GeV 2 and asymptotically remains at larger 
values than for the neutron. 

In Fig. we show our result for f 2 as extracted from Eq. fTOj) . To calculate f 2 , we need 
the full Q 2 dependence of Ti consisting of a sum of elastic and inelastic contributions. The 
elastic contribution is given by Eq. (|2~T]) and evaluated using the most recent experimental 
information on the nucleon form factors as detailed above. The inelastic contribution can 
e.g. be obtained by using the interpolating formula of Eq. (f5T?|) . Although this gives a decent 
description of both proton and neutron data, a better fit for the proton data is obtained by 
using as interpolating formula : 

T ^nel m = T Wres m + ^ as ^ ( , (62 ) 

where is as defined in Eq. (IfiOj) . with I[ es,p (0) = —0.67 obtained from the MAID model. 
Furthermore, we fix the asympto tic v alue by using )p = 0.118 at Q 2 = 5 GeV 2 , as obtained 
from the global analysis of Ref. 3|. This value then fixes the asymptotic coefficient r^ as 
in the interpolating formula for Y^ me \ as well as the flavor singlet axial charge in the 

17 



twist-2 part T ( ±^ tw _ 2 of Eq. (j3J). The shape parameter n in Eq. is taken as jj, ~ 0.74 GeV, 
which gives a very good description of all proton data for 1^ inel as seen from Fig. [7| Note 
that we merely consider the functional form of Eq. (jo^|) as a convenient representation of 
the inelastic data. 

For the neutron, we use the interpolating formula of Eq. ()59j) . with J[ es ' n (0) = —0.48 



obtained from the MAID model. Furthermore, we use as asymptotic value T\ 



(1) as,n 



-0.038, 



and mass scale u ^ 0.57 GeV. Note that in order to get a good description of the JLab/Hall 
A data of Ref . |37J , we use an asymptotic value which is somewhat smaller in absolute value 
than the value rj^ n (Q 2 = 5 GeV 2 ) = -0.058 ± 0.005 ± 0.008 as obtained from the global 



f. Q. 



analysis of Ref. 

Furthermore in Fig.0 we also show separately the twist-2 part T ( [\ w _ 2 , and the term (0.2 + 
4d 2 ) -Mff/(9Q 2 ), which is seen to be negligibly small for both proton and neutron. Therefore 
f 2 is dominated by the difference T<i )el + Y { l )mel - T ( ^ tw _ 2 . The quantity f 2 ■ 4M 2 r /(9Q 2 ) 
evaluated using Eq. (jlOJ) is shown by the thick solid curves in Fig. [7| We can then extract 
f 2 from the linear regime in the 1/Q 2 plot. Because the slight non-linear structure observed 
in the thick solid curves in Fig. is within the error bars of the data, we make a linear fit 
by including all data for Q 2 in the range 0.5 — 2 GeV 2 , in particular the new JLab data 
for both proton and neutron. This is shown by the shaded bands in the same figure, which 
fully accomodate the thick solid curves in the mentioned Q 2 range. From the slopes of these 
bands, we can then extract f 2 for Q 2 in the range 0.5 - 2 GeV 2 as : 



/I 



0.15 0.18, 
-0.026 -> -0.013. 



(63) 



The twist-4 matrix elements f 2 have been extracted before in Ref. ll| from the data of 
Ref. P at Q 2 = 1 GeV 2 , resulting in the values 51 1 : 



tp 

J2 
J2 



0.10 ±0.05 
0.07 ±0.08 



(64) 



Comparing our results of Eq. (|6l3|) with those of Ref. (llj, we see that they are compatible for 
the proton but that we extract a small negative value for the neutron, whereas in Ref. ^j] 
a positive result was extracted. However, given the uncertainty of the neutron data both 
values of f 2 might well be consistent with zero. One sees from Fig. [7[ that for the neutron 
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there is a partial cancellation between the elastic and inelastic contribution to , and the 
result for f 2 is practically completely given by the resonance contribution Y^ res , which is 
estimated here in the MAID model and gives a good description of the available neutron 
data as we have seen before. This mightpartly explain why we extract a negative value for 
f 2 , compared with the analysis of Ref. [111 ]. 

We can compare these phenomenological values for f 2 with several model estimates. In 
Ref. [3], an estimate was given within QCD sum rules which yielded as values : /f = 



2l| . the instanton vacuum picture was 



-0.037 ± 0.006 and / 2 " = -0.013 ± 0.006. In Ref. 
used to estimate matrix elements of quark gluon operators. In particular it has been shown 
in this approach that the matrix element d 2 is suppressed relative to f 2 as : 

d2 rP ^«l, (65) 



/2 \R, 

where p/R ~ 1/3 is the average instanton density. Using the ST/ (3) symmetric case, the 



twist-4 matrix elements f 2 were obtained in Ref. |2l| as : f 2 = —0.046 and f 2 = +0.038. We 
notice that in the approach of Ref. j^], proton and neutron values are of comparable size 
and opposite sign. In both the QCD sum rules and instanton vacuum models, the proton 
value is smaller and of opposite sign compared with the phenomenologically extracted value 
of Ref. ^j] and with the value Eq. ([63)1 of this work. 

It has been suggested that the twist-3 moment d 2 and the twist-four moment f 2 are related 
to the response of the color electric (xe) and magnetic (xb) fields to the polarization of the 
nucleon in its rest frame j3, 4t|, defined as : 

(PS\^gBil>\PS) = X b2M 2 n S, 
(PS\^ ax gEip\PS) = X e2M 2 n S, (66) 

where P is the nucleon momentum and 5* the projection of its spin vector S. Furthermore 
E (B) are the color electric (magnetic) fields respectively, and g is the strong coupling 
constant. Furthermore, the moments d 2 and f 2 can be expressed in terms of the gluon-field 
polarizabilities as : 

1 . 

d 2 = - (Xe + 2xb) , 

/2 = (Xe - Xb) ■ (67) 

Since the experimental value of d 2 is of order 10~ 3 , we see from Eq. f[63[) that the predicted 
central values of f 2 are larger than those of d 2 by about a factor 50 for the proton. These 
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findings agree, at least qualitativel y, w ith estimates using QCD sum rules [23j] and based 
on the instanton vacuum approach j^], but less so with the predictions of bag models jl^ . 
Within the large uncertainties of all existing predictions, we may therefore conclude that 
^2^/2- This observation can be combined with Eq. (J67|) to yield 

2 

XE ~ +g fl , 

Xb « -\h ■ (68) 

In particular a positive value of /2, such as found in both the phenomenological extraction of 
Eq. and Eq. (jfHJ) for the proton, leads to a negative value of xb, i-e., color diamagnetism. 



VI. CONCLUSIONS 



In this work we have studied higher moments of nucleon spin structure functions in view 
of recent data at intermediate values of Q 2 . In particular, we evaluated the generalizations 
to arbitrary momentum transfer Q 2 of the third moment d 2 of the twist-3 part of the nucleon 
spin structure function g2, as well as the third moments and I2 of g\ and gi respectively. 
We have shown that the physical interpretation of d,2, Ii , and I2 & t arbitrary values of Q 2 
is given through nucleon generalized (i.e. Q 2 dependent) spin polarizabilities. These higher 
moments in the Bjorken variable x contain appreciable contributions from the resonance 
region. We therefore evaluated these moments in heavy baryon chiral perturbation theory 
(HBChPT) at order 0(p A ) and in a unitary isobar model (MAID). 

The ChPT results were found to rise strongly with Q 2 , and display a Q 6 dependence at 
low Q 2 values, proportional to the forward spin polarizabilities at the real photon point. 
For d,2, the 0(p 3 ) HBChPT results are in good agreement with the phenomenological MAID 
estimate up to Q 2 ~ 0.25 GeV 2 . However there is a large difference between the 0(p 3 ) 
and 0(p A ) results. This difference originates from the known large difference between the 
0(p 3 ) and 0{p i ) HBChPT results for the forward spin polarizability 70, for which the chiral 
expansion is poorly converging. 

The phenomenological unitary isobar prediction using the MAID model shows an ex- 
cellent agreement with recent neutron data from JLab/Hall A at intermediate Q 2 values. 
The good description is found in part to be due to a sizeable contribution from the ctlt 
photoabsorption cross section. Furthermore, we performed a multipole expansion of these 
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higher moments to gain some additional insight in the dominant absorption mechanisms. 

(3) 

We found that both d 2 and are dominated by an interplay of s- and p-wave multipoles 
up to Q 2 PS 4 GeV 2 , and observed a striking difference between the proton and neutron, 
which is essentially due to the excitation of the S f n(1535) resonance. The helicity amplitudes 
for the photoexcitation of the Su on the proton are about a factor 4 larger as compared to 
the neutron, and this resonance has a very hard form factor. This leads to the result that 
both d 2 and T^ p change sign around Q 2 m 1.25 GeV 2 , and asymptotically remain at much 
larger values than in the case of the neutron. It will be interesting to test this prediction by 

(3) 

extracting d 2 and 1^ from proton data, and to see if one obtains a similar good description 
as for the neutron. 

Besides the twist-3 matrix element d 2 , the 1/Q 2 suppressed term in the twist expansion 
of the first moment of g\ also contains the twist-4 matrix element fa. By including 
all data for for Q 2 in the range 0.5 — 2 GeV 2 , in particular the new JLab data for 
both proton and neutron, we extracted fa in this Q 2 range as : f 2 : 0.15 — > 0.18, and 
/ 2 n : -0.026 -> -0.013. 

The values of d 2 and fa enter in the response of the color electric and magnetic fields to 
the polarization of the nucleon in its rest frame. Therefore our numerical estimates yield 
phenomenological predictions for gluon field polarizabilities, which are new nucleon struc- 
ture information. 

In summary, the sum rules which we studied in this work as function of Q 2 from the 
real photon point to the region of deep inelastic scattering, are very interesting observables 
to interpolate between the non-perturbative and perturbative regimes of QCD. As these 
higher moments do not require large extrapolations into unmeasured regions, they are ideal 
observables to be measured at intermediate momentum transfers. 
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Figure 1: Q 2 dependence of the moment di for proton (upper panel) and neutron (lower panel). 
The dashed-dotted curve is the elastic contribution to di according to Eq. The other curves 

represent the inelastic contributions to d,2- Solid curves : MAID estimate for the tt channel; 
dotted curves : 0(p 3 ) HBChPT; thick (upper) dashed curves : C(p 3 ) + 0(p A ) HBChPT; thin 
(lower) dashed curves : 0(p 3 ) HBChPT with 0(e 3 ) A contribution added. The JLab/Hall A data 
(diamonds) are from Ref. J4| (inner error bars are statistical errors only, outer error bars include 
systematical errors). The SLAG data (circles at = 5 GeV*) are from Ref. Q. 



24 



0.01 
0.005 


-0.005 























7t + TOT + Tl 








71 


r (3)p 
1 1 



0.5 



0.015 
0.0125 
0.01 
0.0075 
0.005 
0.0025 


-0.0025 
-0.005 
-0.0075 



0.5 



1.5 2 2.5 3 3.5 4 













1 7 









1.5 2 2.5 3 3.5 4 

Q 2 (GeV 2 ) 



Figure 2: Q 2 dependence of the moment for proton (upper panel) and neutron (lower panel). 
The dashed-dotted curve is the elastic contribution to 1^ according to Eq. (|21j) . The other curves 

(3) 

represent the inelastic contributions to T\ ' . Thin solid curves : MAID estimates for the 7r and 
7r + 7T7T + rj channels (for proton) as indicated on curves; upper thick solid curve (for proton) 
: total estimate according to Eq. (|59|): dotted curves : 0(p 3 ) HBChPT; thick dashed curves : 
0{p 3 ) + 0(p 4 ) HBChPT; thin dashed curves : 0(p 3 ) HBChPT with 0(e 3 ) A contribution added. 
The upper shaded band is the evaluation using the DIS structure function g\ as extracted from 
experiment, according to the analyis of Ref. [32|. The lower shaded band is the corresponding 
analysis for the contribution from the region W < 2 GeV (i.e. resonance region) to the DIS 
parametrization. The size of the bands represents the corresponding (la) error estimates as given 



by Ref. 



32] 
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Figure 3: Q 2 dependence of the moment T^ for proton (upper panel) and neutron (lower panel). 
The dashed-dotted curves are the elastic contributions to T^' according to Eq. (|24|). The other 

(3) 

curves represent the inelastic contributions to rj> • Solid curves : MAID estimates for the tt 
channel; dotted curves : 0{p 3 ) HBChPT; thick dashed curves : 0(p 3 ) + 0{p A ) HBChPT; thin 
dashed curves : 0(p 3 ) HBChPT with C(e 3 ) A contribution added. 
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Fi gure 4: Different photon absorption cross section contributions to d,2 (upper panel) and 
(lower panel) for the neutron. The solid curves are the total MAID estimate for the tt channel. 
The dashed (dotted) curves are the contributions from &tt (plt) separately in Eqs. ((181) and (|54|). 
Data for cfo as in Fig. ^ 
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Figure 5: Multipole content of the moment di for proton (upper panel) and neutron (lower panel). 
The solid curves are the total MAID estimate for the tt channel (same as solid curves in Fig. 
The other curves represent different partial wave contributions. Dashed-dotted curves : results for 
the sum of s- and p- waves; dashed curves : results for only s- waves; dotted curves : results for 
s-waves without 5n(1535) resonance contribution. Data as in Fig. ^ 
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Figure 6: Multipole content of the moment T\ for proton (upper panel) and neutron (lower panel). 
The solid curves are the total MAID estimates for the tt channel (same as corresponding solid curves 
in Fig. [2J). The other curves represent different partial wave contributions. Dashed-dotted curves 
: results for the sum of s- and p-waves; dashed curves : results for only s-waves; dotted curves : 
results for s-waves without 5n(1535) resonance contribution. 
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Figure 7: Q 2 dependence of fa ■ / (9Q 2 ) (thick solid curves in the shaded bands) for proton 
(upper panel) and neutron (lower panel) as extracted from Eq. (|1U|) . We show separately the 



elastic T^ el (thin dashed-dotted curves), total inelastic Y\ 1>inei (thick dashed-dotted curves), and 
resonance ( W < 2 GeV ) (thin solid curves) contributions to We also show the twist- 2 

part r^ U) _ 2 (dashed curves) and (02 + 4c?2) ■ Afjy / (9Q 2 ) (dotted curves) which enter on the rhs 
of Eq. ((31) for T± . The shaded bands are a linear fit to extract fo as described in the text. The 



(1) inel 



proton data for T 



(1) inel 
1 



are from JLab/CLAS |2i (circles), SLAC 



4jJ (triangles). The neutron data for pW* nei are from JLab/HallA 
at Q 2 = 5 GeV 2 (stars) correspond with the global analysis of Ref. 
statistical errors only, outer error bars include systematical errors. 



(diamonds), and HERMES 
(circles). The data points 



451 ] . Inner error bars are 
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